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Abstract 

In the Nagel-Schreckenberg model of vehicular traffic on single-lane 
highways vehicles are modelled as particles which hop forward from one 
site to another on a one dimensional lattice and the inter-particle interac- 
tions mimic the manner in which the real vehicles influence each other's 
motion. In this model the number of empty lattice sites in front of a parti- 
cle is taken to be a measure of the corresponding distance-headway{DYl) . 
The time-headway(TH.) is defined as the time interval between the depar- 
tures (or arrivals) of two successive particles recorded by a detector placed 
at a fixed position on the model highway. We investigate the effects of 
spatial inhomogeneities of the highway (static hindrances) on the DH and 
TH distributions in the steady-state of this model. 



PACS. 05.40.+j - Fluctuation phenomena, random processes, and Brownian 
motion. 

PACS. 05.60.+W - Transport processes: theory. 
PACS. 89.40. +k - Transportation. 
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1 Introduction: 

The continuum models of traffic flow [f -3] are analogues of the " hydrodynamic" 
models of fluid flow while the kinetic theories of vehicular traffic [4-6], which 
are extensions of the kinetic theory of gases, and the car-following models [7-10] 
as well as the discrete "particle-hopping" models [11-17] are analogues of the 
"microscopic" models of interacting particles commonly studied in statistical 
mechanics. In this paper we focus our attention on a specific particle-hopping 
model, namely, the Nagel-Schreckenberg (NS) model [11] of vehicular traffic 
on idealized single-lane highways; this model may be regarded as a model of 
interacting particles driven far from equilibrium and the dynamical phenomena 
exhibited by this model of traffic may be treated as problems of non-equilibrium 
statistical mechanics. 

The distance-headway(DH) is defined as the distance from a selected point 
on the lead vehicle (LV) to the same point on the following vehicle (FV). 
Usually, the front edges or bumpers are selected [18]. The time-headway(TH.) is 
defined as the time interval between the departures (or arrivals) of two successive 
vehicles recorded by a detector placed at a fixed position on the highway [18]. 
The DH and TH are not merely of academic interest to statistical physicists, but 
are also of practical interest in traffic engineering where these two distributions 
are regarded as important characteristic of traffic flow [18,19]. For example, 
larger headways provide greater margins of safety whereas higher capacities of 
the highway require smaller headways. 
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Recent investigations by several groups of statistical physicists [20-24] have 
helped in gaining insight into the nature of the "jamming transition" (i.e., the 
transition from the " free- flowing" dynamical phase to the "jammed" dynamical 
phase) in the NS model. On the other hand, simultaneously NS model is 
being generalized and extended [25-37] to capture more and more details of 
vehicular traffic flow so that the generalized/extended model may, ultimately, 
find practical use in traffic engineering [38-39] . Extending our earlier numerical 
works on the DH distribution in the steady-state of the NS model, we point out 
here interesting similarities as well as crucial differences between the "jamming 
transition" and the gas-liquid phase transition in a simple fluid in equilibrium. 
Moreover, in this paper we also present a simple analytical calculation of the 
TH distribution in the steady-state of a special case of the NS model using an 
interesting quantity introduced in ref.[15]. Furthermore, we study the effects of 
spatial inhomogeneities of the highway (static hindrances) on the DH and TH 
distributions in the steady-state of the NS model. 

The NS model and some of its most important features are summarized in 
section 2. We report numerical results on the DH distribution in the steady- 
state of the NS model, in the absence as well as in the presence of hindrances, in 
section 3 where we also draw attention to the analogies and differences between 
the jamming transition in the NS model and the gas-liquid transition. Our 
calculations of the TH distribution in the steady-state of the NS model are 
given in the section 4, where we also study the effects of hindrances on this 



3 



distribution. Finally, we summarize our main results and conclusions in section 
5. 

2 The Models: 

In the NS model a lane is represented by a one-dimensional lattice of L sites. 
Each of the lattice sites can be either empty or occupied by at most one "ve- 
hicle" . If periodic boundary condition is imposed, the density c of the vehicles 
is N/L where N(< L) is the total number of vehicles. In the NS model [11] 
the speed V of each vehicle can take one of the V max + 1 allowed integer values 
V = 0, 1, Vmax- Suppose, V n is the speed of the n-th vehicle at time t. At 
each discrete time step t — ► t + 1, the arrangement of N vehicles is updated in 
parallel according to the following " rules" : 

Step 1: Acceleration. If, V n < V max , the speed of the n-th vehicle is increased 
by one, i.e., V n ->V n + l. 

Step 2: Deceleration (due to other vehicles). If d n is the gap in between the 
n-th vehicle and the vehicle in front of it, and if d n <V n , the speed of the n-th 
vehicle is reduced to d n — 1, i.e., V n — > d n — 1. 

Step 3: Randomization. If V n > 0, the speed of the n-th vehicle is decreased 
randomly by unity (i.e., V n — > V n — 1) with probability p (0 < p < 1); p, the 
random deceleration probability, is identical for all the vehicles and does not 
change during the updating. 
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Step 4-' Vehicle movement. Each vehicle is moved forward so that X n — > X n + V n 
where X n denotes the position of the n-th vehicle at time t. 

The specific update rule of the NS model requires a nonvanishing braking 
probability p for the model to yield a realistic description of traffic flow [12] and, 
thus, the NS model may be regarded as stochastic cellular automata [40]. Effec- 
tively free flow of traffic takes place when the density of vehicles is sufficiently 
low whereas high density leads to congestion and traffic jams. One of our aims 
is to point out some similarities as well as differences between the "jamming 
transition" in the NS model [20-24], which is a non-equilibrium driven system, 
and the gas-liquid phase transition in a fluid in equilibrium. 

The asymmetric simple exclusion process (ASEP) is one of the simplest mod- 
els of driven systems of interacting particles; it is often treated as a carricature 
of traffic flow. The relations between this model and the V max = 1 limit of 
the NS model has been elucidated in the literature [24]. Some dramatic effects 
of quenched disorder (static hindrances) on the steady state of this model have 
been investigated [41,42]. Our main aim in this paper is to investigate the effects 
of static hindrances on the DH and TH distributions in the NS model. 

Emmerich and Rank [32] introduced an extra step of update rule before 
all the other steps of updating in the NS model to mimic the effects of static 
traffic hindrance: a segment of length Lhind on the highway is identified as the 
hindrance and the speed of all the vehicles found within that segment of the 
highway are reduced to half of their current speed. A more general model of 
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traffic flow in the presence of " quenched disorder" in the highway was formulated 
by Csahok and Vicsek [33] in terms of an " inverse permeability" associated with 
each site. We shall use the simpler rule introduced by Emmerich and Rank [32]. 

For the convenience of our analytical calculations, following Schrcckenberg 
et al.[14], we assume the sequence of steps 2 — 3 — 4 — 1, instead of 1 — 2 — 3 — 4; 
the advantage is that there is no vehicle with V = immediately after the 
acceleration step. Consequently, if V max = 1, we can then use a binary site 
variable a to describe the state of each site; a = represents an empty site and 
cr = l represents a site occupied by a vehicle whose speed is unity. 

3 Distance- Headway Distributions and the Na- 
ture of the Jamming Transition 

Using the sequence of steps 2 — 3 — 4 — 1, as explained earlier, an n-cluster con- 
figuration in the steady-state of the NS model is represented by (cti, 02, <r„). 
The number of empty lattice sites, n, in front of a vehicle is taken to be a 
measure of the corresponding DH. Within the 2-cluster approximation, the DH 
distribution, T"^{n), in the steady-state of the NS model with V max = 1 has 
been calculated analytically following the methods of the site-oriented mean- 
field (SOMF) theory [14]. In this approximation, 

7>*(0)=C(1|1) (1) 
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and 

^W=C(l|0){C(0|0)r- 1 C(0|l), for n>\ (2) 

where, C gives the 2-cluster steady-state configurational probability for the ar- 
gument configuration and the underlined imply the conditional, as usual. The 
expressions for the various Cs are given by [14,24] 

C(0|0) = C(0|0) = 1 - - d (3) 

C(1|0) = C{Q\1) = V - (4) 
C(Q\l) = C(1|0) = | (5) 

caii) = c(i|i) = i-| (6) 

where 

y = Y q ( 1 - yfi - 4 H - ( 7 ) 

q = l—p and d = 1 — c. It turned out that the 2-cluster approxmation is exact 
for V m ax = 1 [14]. Therefore, in terms of c and p, the exact DH distribution in 
the steady-state of the NS model with V max — 1 is given by 

7^(0) = 1 - (y/c) (8) 

and 

V dh {n) = {y 2 /(cd)}[l - (y/dW- 1 for n > 1 (9) 

where, for the given c and p, y can be obtained from equation (7). The same 
exact DH distribution (8-9) has also been derived independently [15] within the 
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framework of car-oriented mean-field (COMF) theory. The DH distributions 
in the NS model for all V max > 1 have been computed numerically by carrying 
out computer simulation [24]. 

The DH distribution in the steady-state of the NS model simultaneously ex- 
hibits two peaks over an intermediate regime of the vehicle densities, provided 
V m ax > 1 and p is sufficiently large. In contrast, when V max = 1, no such two- 
peak structure is exhibited at any density, for any p, by the DH distribution 
which is given by the exact expression (8-9). We have now estimated, as func- 
tions of p, the densities, c s and eg, corresponding to the smallest and the largest 
densities where two-peak structure is exhibited by the DH distribution in the 
NS model (2 < V max < 5). c s and a are plotted against 1 — p for V max = 2 in 
fig.l. 

The occurrence of the two-peak structure for V max > 1 has been interpreted 
[24,43] as a signature of "two-phase coexistence"; the two coexisting phases, 
namely, the " free- flowing phase" and the "jammed phase" being the analogues of 
the gas phase and the liquid phase, respectively, of a simple fluid in equilibrium. 
In this scenario, for a given p (and given V max ), it is tempting to identify c s as 
the analogue of the density corresponding to the onset of two-phase coexistence 
in a fluid while eg would be the corresponding largest density up to which these 
two phases continue to coexist at a given temperature; accordingly, the regime 
of density in between these two curves would be the analogue of the two-phase 
coexistence region for a fluid and the uppermost tip of this region could be 
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interpreted as an analogue of the critical point. Further, this analogy would 
suggest that 1 — p is the analogue of temperature; increase of p (i.e., decrease 
of 1 — p) drives the system towards " condensation" . 

In spite of this apparent analogy, there are crucial differences between the 
"jamming transition" and the gas-liquid transition in a simple fluid. First of 
all, the fluid is in equilibrium whereas the model vehicular traffic under consid- 
eration is in a non-equlibrium steady-state. More important difference is that 
the size of the jams remain finite even in the limit of infinite system size [24] 
and, therefore, cannot be identified as a true dynamical phase. 

When the DH distribution exhibits two peaks simultaneously, the peak at 
vanishing DH reflects the fact that the total number of vehicles held up simul- 
taneously in the various jams on this model highway is a finite fraction of all 
the vehicles on this highway; these jams arise from spontaneous fluctuations in 
the driven system of interacting particles and, therefore, can appear anywhere 
in the system. In contrast, a hindrance can induce recurring jams at (and near) 
its own location, as demonstrated in fig. 2, because of the bottlenecks against 
traffic flow created by it. The DH distributions for 5 different values of L hind 
are plotted in fig. 3, all for the same density c = 0.05 and for the same p = 0.5. 
In this figure the absence of two-peak structure for Lhind = indicates that, 
for p — 0.5, the density c = 0.05 is not high enough to have a finite fraction 
of the vehicles, simultaneously, held up in various jams created by spontaneous 
fluctuations. However, for the same set of values of c and p, two-peak structure 
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occurs in the DH distribution with the increase of Lhind] this is a consequence 
of jams caused by the bottleneck effect of the hindrance. 

4 Time-Headway Distribution: 

There are several earlier papers, published by statisticians and traffic engineers, 
where the form of the TH distribution has been derived on the basis of heuristic 
arguments [44] . In this section we begin by presenting a simple derivation of this 
distribution for V max = 1 using a quantity introduced in ref.[15] in the context 
of the COMF theory. However, we could compute the corresponding TH dis- 
tributions for V max > 1 [45] and the effects of hindrances on these distributions 
only through computer simulation. 

We label the position of the detector by j = 0, the site immediately in front 
of it by j = 1, and so on. The detector clock resets to t = everytime a vehicle 
leaves the detector site. We begin our analytical calculations for V max — 1 by 
writing V th (t), the probability of a time headway t between a LV and the FV, 
as 

t-i 

v th {t) = Y,P(ti)Q'(t-ti\t 1 ) (io) 
ti=i 

where P{t\) is the probability that there is a time interval t\ between the de- 
parture of the LV and the arrival of the FV at the detector site and Q'(t — 1\ \t\) 
is the conditional probability that the FV halts for t — t\ time steps when it 
arrived at the detector site t\ time steps after the departure of the LV. 
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Suppose, g(t) is the probability that a vehicle moves in the next (i.e., in the 
t+l-th time step). Therefore, g(t) = l—g(t) is the probability that a vehicle does 
not move in the next time step. It has been shown [15] that g = q[l — V dh (0)]. 
Using equation (8) for V dh (0) we get 

g = qy/c (11) 

and, hence, g = 1 — qy/c, where y is given by the equation (7). Moreover, since 
y satisfies the equation qy 2 — y + cd = 0, it follows that 

(l-qy/c)(l-qy/d)=p. (12) 

In order to calculate P{t\) we need to consider all those spatial configurations 

at t = from which the FV can reach the detector site within t\ steps. For all 

configurations with t\ > n, t\ — 1 time steps elapse in crossing n — 1 linkss (as 

the last link is crossed certainly at the last time step). Thus 

*i 

P(h) = J2^W l p tl - n ^CVi (13) 
n=l 

where 

n(n)=C(l|0){C(0|0)}- 1 . (14) 

and, hence, we get 

P(h) = c(i\o) q [cmv + p]* 1 ' 1 = f (i - f ) tl_1 (is) 

which is the exact analytical expression for P{t\). 

Next, we calculate Q(t — t\ \ t\) by expressing it in terms of g and g. When 
the FV arrives at the detector site exactly t\ time steps after the departure of 
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the LV, the LV can be at any of the sites labeled by 1 , • • • , 1 1 + 1 . The probability 
that the LV stays at the site T' is g tx and, therefore, the probability that the 
LV is not at the site T' is 1 - g tx . 

If the LV is not at site T' then the probability that the FV halts at the 
detector site for exactly t — t\ time steps is p t ~ tl ~ 1 q because it should stop due 
to randomisation for exactly t — t\ — 1 steps and move at the last step. Hence, 
the contribution to Q'(t — t\ \ t\) when the LV is not at site T' is 



tn„*-ti-i. 



(I) 



On the other hand, when the LV is at site T' then it will have to move so that 
the FV is able to leave the detector site after t — 1\ steps. Suppose the LV moves 
from T' after k steps (k varies from 1 to t — 1\ — 1). Then the FV will have to 
stay at the detector site for next t — t x — 1 — k steps due to randomisation and 
move in the last step. Hence, the contribution to Q'(t — t\ \t{) when LV is at 
site T' is 



t-ii-l 



g H gq{ E 9 k - 1 P t - tl - 1 - k ) = g tl gq 



k=i 



g-p 



(ii) 



Therefore, combining (I) and (II) we get 



t-ti-ii 



Q\t - h | h) = (1 - g^p^^q + g u gq 



g-p 



(16) 



Finally, substituting the exact expressions (15) and (16) into (10) we get 



V th (t) 



qy 



c-y 



1 1— l 



{i - (qy/c)} 

qy qy 



qy 



{i - (qy/d)} 



t-i 



c-y d-y 



d — y 

p*- 1 - q\t - iy- 2 . (17) 
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where, for the given c andp, y can be obtained from equation (7). The invariance 
of the distribution (17) under the interchange of c and 1 — c, is a manifestations 
of the well-known " particle-hole" symmetry in the problem which breaks down 
for all V max > 1 [46]. 

The flux q of the vehicles can be written as q = N/T where T = U 1S 

the sum of the time headways recorded for all the N vehicles. Therefore, one can 
rewrite q as q = 1/T av where T av = (l/N) J2i U is the average TH. Therefore, 
T av is expected to exhibit a minimum at c = c m with the variation of density 
c of the vehicles. We observed that the trend of variation of the most-probable 
TH, T mp , with c is similar to that of T av with c [45]. Moreover, because of the 
particle-hole symmetry, the T mp versus c curve is symmetric about c = 1/2 in 
the NS model with V max — 1. But this symmetry is lost when V max > 1. The 
fact that T mp (and T av ) versus c curve exhibits a minimum is consistent with 
one's intuitive expectation that both at very low and very high densities there 
are long time gaps in between the departures of two successive vehicles from a 
given site. 

Two typical sets of curves showing the effects of the hindrances on the TH 
distribution in the steady-state of the NS model are shown in fig. 4. At low 
density of the vehicles, increase of the length of the hindrance leads to significant 
broadening of the distribution although its effect on the magnitute of the most 
probable TH is weak (see fig. 4a). On the other hand, at moderate and high 
densities, the broadening of the TH distribution caused by increase of L hind is 
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much weaker (see fig. 4b). 



5 Summary and Conclusion: 

In this paper we have investigated two important characteristics of traffic flow 
on highways, namely, DH and TH distributions, starting from models that in- 
corporate explicitly vehicle-vehicle and road-vehicle interactions. We have com- 
pared and contrasted the effects of static hindrance on the highway and those 
of the LV, which may be viewed as a dynamic hindrance, against the forward 
movement of the FV. In the absence of any static hindrance jams can appear 
anywhere in the system because of spontaneous fluctuations and the dynamic 
hindrance caused by a slower LV on the FV. On the other hand, in the presence 
of a static hindrance, jams take place at (and near) the hindrance, in addition 
to those formed by spontaneous fluctuations, because of bottleneck it creates 
against the traffic flow. 

We have extracted some informations on the " structures" in the spatial or- 
aganization of the vehicles on the highway from the DH distribution. We have 
demonstrated the effects of the bottleneck against traffic flow created by static 
hindrances on the DH distribution thereby elucidating the physical meaning 
of the two-peak structures observed earlier in the DH distribution over some 
intermediate regime of density when V max > 1 and p is sufficiently large. 

The TH distributions in the steady state of the NS model and the trend 
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of their variation with density are in good qualitative agreement with the cor- 
responding empirical data [18]. Our results demonstrate that, in spite of its 
simplicity, the NS model captures the essential qualitative features of the TH 
distribution of vehicular traffic on highways. Besides, the extent of broadening 
of the TH distribution caused by a hindrance depends on the density of the 
vehicles; the higher is the density the larger is the broadening. Hoever, at all 
densities, the effect of the hindrance on the magnitude of the T mp is very weak. 

In this paper we have considered only the original version of the NS model 
for vehicular traffic on single-lane highways [11] and an extended version that 
includes a specific type of inhomogeneities of the highway (namely, static hin- 
drances) [32]. We have not attempted any direct quantitative comparison of 
our results with the corresponding empirical data from highway traffic because, 
we believe, such comparisons will be possible only after several of the realistic 
generalizations and extensions [25-37], proposed recently in the literature, are 
incorporated in the model. Results of our ongoing works in this direction will 
be published elsewhere [47]. 
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Figure Captions: 

Fig.l: ci and c s in the NS model are plotted against 1 — p when V max = 2 (see 
the text for the definitions of q and c s ). 

Fig. 2: The "space-time diagram" showing the time-evolution of the traffic in 
the NS model (V max = 5,p = 0.5) in the presence of ten hindrances, each of 
length 5, put randomly along the highway of length 1000. Each of the black 
dots represents a vehicle. 

Fig. 3: The distance-headway distributions corresponding to the vehicle density 
c = 0.05 in the NS model with V max = 5 (p = 0.5) for five different values of 
the hindrance length are compared with that in the absence of hindrance. The 
discrete data points obtained from computer simulation correspond to Lhind = 
0(+),L hind = l(x),L hind = 2(*),L hind = 3(a),L hind = 4(m),L hind = 5(o) 
while the continuous curves are merely guides to the eye. 

Fig. 4: The time-headway distributions in the NS model (V max = 5), in the 
presence of hindrance, for vehicle densities (a) c = 0.1 and (b) c = 0.5. The 
symbols corresponding to Lhind — 0,1,2,3,4 and 5 are identical to those in 
fig.3. 
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